February 2, 2008 16:12 WSPC/INSTRUCTION FILE flavorphysics-07 



International Journal of Modern Physics A 
© World Scientific Publishing Company 



Semileptonic B(Bs, Be) decays in the light-cone QCD sum rules 



Tao Huang* 

Institute of High Energy Physics, Chinese Academy of Sciences, Beijing 100049, China 

Zuo-Hong Li 

Department of Physics, Yantai University, Yantai 264005, China 
Xing-Gang Wu 

Department of Physics, Chongqing University, Chongqing 4 00044i China 

Fen Zuo 

Department of Modern Physics, University of Science and Technology of China, Hefei, Anhui 

230026, China 

Institute of High Energy Physics, Chinese Academy of Sciences, Beijing 100049, China 



Received Day Month Year 
Revised Day Month Year 

Semileptonic B{Bs, Be) decays are investigated systematically in the light-cone QCD 
sum rules. Special emphasis is put on the LCSR calculation on weak form factors with 
an adequate chiral current correlator, which turns out to be particularly effective to 
control the pollution by higher twist components of spectator mesons. The result for 
each channel depends on the distribution amplitude of the the producing meson. The 
leading twist distribution amplitudes of the related heavy mesons and charmonium are 
worked out by a model approach in the reasonable way. A practical scenario is suggested 
to understand the behavior of weak form factors in the whole kinematically accessible 
ranges. The decay widths and branching ratios are estimated for several B{Bc) decay 
modes of current interest. 
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1. Introduction 

Semileptonic B{Bs, Be) decays provide an important gromid to understand and test 
the standard model(SM), and perhaps, a window into new physic beyond the SM. 
To test the SM with the experimental data demands that we have the ability to 
precisely compute the physical amplitudes in QCD theory. In this talk we focus on a 

*Email: huangtao@mail.ihep.ac.cn 

1 



2, 2008 16:12 WSPC/INSTRUCTION FILE flavorphysics-07 



2 Tao Huang, Zuo-Hong Li, Xing-Gang Wu and Fen Zuo 



discussion about how to compute the related hadronic matrix elements in the QCD 
light cone sum rules (LCSR's) which have proved to be a powerful tool to derive 
the desired form factors with a minimal number of phenomenological assumptions. 
Instead of the procedure followed by conventional sum rule calculations, a LCSR 
computation starts with a vacuum-meson correlator, with the related light meson 
being on its mass shell, and operator product expansion (OPE) of the correlator 
is performed in term of the vacuum-meson matrix elements of certain nonlocal 
operators near light cone = 0. Then one parameterizes these nonlocal matrix 
elements in terms of light-cone wavefunctions of the light meson with increasing 
twist. In contrast to the case of the 3-point sum rules (3PSR's), the effective region 
of the momentum transfer can be fixed without any extrapolation. However, it 
is our concern how to effectively control the nonperturbative dynamics embedded 
in the light-cone wavefunctions for enhancing the predictive power of the LCSR 
method. As shown in Ref. , the contributions would, in general, be important 
from some of the higher twist distribution amplitudes, especially twist-3 ones, in 
many situations. A pragmatic prescription 1^ has been put forward by us to reduce 
the contamination by the higher- twists. In the improved LCSR approach, we choose 
to use a certain proper chiral current operator as interpolating fields in the related 
correlators, which make twist-3 wavefunctions disappear in the sum rule results. 
Accordingly, this can enhance considerably one's confidence in applying LCSR's to 
calculate the nonperturbative quantities. 

In this paper, we would like to give the complete discussion of the semileptonic 
form factors for B{Bs,Bc) meson with the chiral current correlators. A model for 
leading twist distribution amplitudes is formulated in the reasonable way, for the re- 
lated heavy mesons and charmonium. Also, a practical scenario is suggested to have 
an all-around understanding of the form factors in the whole kincmatical ranges. 



2. B, Bs — > -P(^) transition form factors 

The hadronic matrix elements for B, Bs PiV) transitions can be parameterized 
in term of the form factors in the following way: 

<Pip)\q7^.b\B^s)ip + q)> = f+{q^)i2p + q)^ + f-iq^)q^, (1) 
<V{p,v)\qi^{l - 75)6|B(.)(p + q)>= - t7j;{mB+mv)A,{q') + i{rj*q)^±^^^^^±^ 

+ 1) \ : , (2) 

iJiB + mv TUB + mv 

where q and 77^ are the momentum transfer and polarization vector of the vector 
meson, respectively. 

For performing a LCSR computation of the form factors for B{Bs) P, we 
choose the following correlator 11^ (p, q) with the chiral current 131^ 

n^{p,q)^l I d^xe'^^<P{p)\T{q^{x)-f^{l+-f^)b{x)~b{m^+lMm^> (3) 
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By contracting the 6-quark operators to a free propagator, we get leading twist 
contribution to the correlator, 

7^^^^'"'^''^^3^<^(^^)l^9"i(^)^^^5«2(o)|o> 

= 2fpmi,Pfj, / du — 2 7 + higher twist terms (4) 

where we have substituted the definition of the leading twist distribution amplitude 
ipp{u): 

<P{p)\Tqi{x)jf,j5q2{0)\0> = -ip^fp f due"'P-^(^p(w)+ higher twist terms. (5) 

Jo 

Because of the chiral nature of the selected correlator ([3]), only the leading non- 
local matrix element <P(p)|gi(x)7^75g2(0)|0> contributes to the correlator, while 
those with leading twist-3 <P{p)\qi{x)ij5q2{0)\0> and <P(p)|gi(a;)cr^,y75q2(0)|0> 
disappear from the sum rule. Including the effect of the background gluon field 
by writing down a full 6-quark propagator, it is shown that all the twist-3 terms 
actually vanish in the OPE of this special correlator. Twist-4 and yet higher twist 
terms are found to be numerically small, and the resulting uncertainties can always 
be neglected comparing with those of the approach itself. Matching Eq. Q to the 
hadronic representation of the correlator and omitting QCD radiative corrections, 
we get the following simple LCSR result for f+{q^) 

^ ^e'"-/^^^ du^ exp [- '"'""iCr"^^ ] + twist-4 terms 

(6) 



with u = 1 — u and 
Ap 



where denotes the corresponding Borel variable, and sq the threshold parameter. 
When setting Mp — 0, we get a simpler Ap = "p . As a byproduct, we also find 

^0 1 

the relation f-{q^) = —f+{q'^), which holds up to twist-4 terms. We would like to 
stress that the resulting form factors are just valid for < < m^ — 2rnbAQCD- 

Form factor f^{q^), as an illustrative example, has been calculated in terms 
of the above expression in the cases of B ^ TT and -> X El. QCD radiative 
corrections have been added to the twist-2 part in the B n case with a neg- 
ligibly small numerical impact. By applying LCSR's to do such calculations, what 
is meant is that one is assuming soft exchanges to predominate in this transition. 
Thus, the semileptonic from factors in the whole kinematical range can be deter- 
mined by combing the LCSR results with those of lattice QCD, which turn out to 
be available near the smallest recoil. While an alternative choice is made, which 
assumes hard exchanges to dominate at larger recoil for which the PQCD approach 
is applicable, a combined use should be in order of the three different approaches!^, 
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to fully understand the behavior of the form factor in the whole kinematic regions. 
A similar discussion on the case oi B —i- K form factor is presented in Ref.EI. 

It is straightforward to extend our discussion to the B ^ D case. In Ref.l^, 
the B ^ D form factor is calculated using three different models for the D meson 
distribution amplitude. In order to compare with the results of heavy quark effective 
theory (HQET), the form factor has been redefined according to 



Ml') = 



niB + mo 
2y/mB mo 



Tb^d{v ■ v' 



(8) 



The numerical results, which correspond to the region < < (ms — moY , are 
shown in Fig[T] At zero recoil JC-iC^(l) = 1.02 (using model III in Ref. '^), which is 
in a good agreement with the evaluation obtained using the heavy quark symmetry: 
^b^d{^) — 0.98±0.07 ^. In the larger recoil region 1.35 < v-v' < 1.59, the yielded 
results are consistent with those of pQCD^Sl. Therefore, a full understanding of the 
dynamics involved in B ^ D transition may be obtained by combining the three 
different approaches — HQET, LCSR's and pQCD, which are complementary to 
each other. 

For B{Bs) V modes, we choose the following correlator as our starting point: 



n^(p,(?) = -ij d4a;e''--<y(p,77)|r{gi(.T)74l-75)6(x),6i(0)(l + 75)g2(0)}|0> 
= T^T^l - r+(77*g)(2p + q)^ - T-{i^*q)q^ + zr^£^„^^77*"gV. (9) 



form factor 




1.0 1.1 1.2 1.3 



1.5 1.6 



Fig. 1. Tb-^d as a function of the velocity transfer. The thin lines expresses the experiment fits 
with the HQET result at zero momentum transfer as input, where the solid line represents the 
central values, the dashed (dash-dotted) lines give the bounds from the linear (quadratic) fits. The 
thick lines correspond to our results, with the solid, dashed and dash-dotted lines for model III, II 
and I respectively. 



2, 2008 16:12 WSPC/INSTRUCTION FILE flavorphysics-07 



Semileptonic B (Bs , Be) decays in the light-cone QCD sum rules 5 



Also we take the standard definition of the twist-2 and twist-3 distribution amph- 
tudes of the vector meson: 



niv 



Actually there are also two other twist-3 distribution ampHtudeslI^, h^^l'\u), h''^^' (u), 
which are related to the square of light meson mass, and therefore can be neglected 
in comparison with those written above. We also neglected twist-4 and higher terms, 
and the whole three-particle contributions. Following the standard procedure, one 
obtains: 



^1(9') 



fBm?g{mB + mv) 



B 

1 du 
— exp 



^6 — u{q^ — urriy) 



1^1 — q + W^niy 



^^2^ ^ f^mbjniB + my) ^^2^^j^f2 



fsmj: 
du 



— exp 



ml — u{q^ — umy) 



M^), (11) 



(12) 



-A+{q') 



(13) 
(14) 



with 



Ay - [^{sX - 92 - ml^y + Amliml - q^) - (^ - - ml)]/ {2ml), (15) 

The same technic can be used to discuss the rare decay B V^. The electromag- 
netic penguin operator dominates in this transition, and relevant decay amplitude 
reads 



A{B ^ Vj) = Cm^e^^ < V{p,r,)\qa^,{l + j5)q''b\B{p + q) >, 



(16) 



where and q are the emitted photon polarization vector and momentum, respec- 
tively. The constant C depends on the product of CKM matrix elements VtbV^j^ 
and the corresponding Wilson coefficient C7. The hadronic matrix element may be 
parameterized in terms of the form factor T{q^ — 0), 

< V{p, r;)|ga^,(l + j5)q''b\B{p + q) >= 2(ie^,„077*''gV + P ■ qvl - q ■ V*P^.)T{0). 

(17) 

A correlator, which is suitable for our purpose, is 

F^,{p,q) =i J d^xe^'^-- <V{p,7])\Tq-i{x)a^,{l + j5)q''b{x),b{m^ + 75)q2(0)\0> 
= [2ze^.„/3r7*"9V + 2p • qv; ~ 2q ■ r,*p^] F [{p + q)^] . (18) 
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According to the definition in Eq. (|10p we can parameterize the nonlocal ma- 
trix element < V{p,'q)\qi{x)<7^ijq^ {1 +75)g2(0)|0) > using the the leading twist-2 
wavefunction (j}^{u) as the following, 

< ^(p,'7)k~i(a;)o-pi.'7''(l + 75)92(0)10 > =i[{q- V*)Pp - P ■ 1% - i^t^^^ap^f q"p'^] fl 

1 

X / due™f^</._L(w). (19) 



Then we get the final LCSR for the form factor r(0) 

Tiq^ = 0) = 4#e™^/^'^^- exp [^ ^g + z^jlO 

where Ay — Av|q2^Q. Notice that 

T{q'-0)^ 7 (21) 
(m_B + mv) 

With th is sum rule, both B K*j and B (p, lu)j have been discussed in 
Ref. I ^^ l '*^^ ! In j^gf ^jjg fQpjji factor and branching ratio for B (p, ui)^ are 
achieved using two different models for the p meson distribution amplitude: 

yB^P± (0) = 0.335 ± 0.050, Br{B -> ^=^7) = (2.71 ± 1.00) x 10"^ 

for the model introduced by Ball and Braun in Ref. ^31^ and 

rj^B^p^ (0) = 0.272 ± 0.029, Br{B p^-f) = (1.79 ± 0.61) x 10"^ 

for the model presented by Bakulev and Mikhailo in Ref. Comparing the LCSR 
predictions with the recent experimental result 

Br{B ^ p±7) = (1.32t°:|t(stat)t°:i°(syst)) x 10-^ 

we can determine which model can better describe the p meson. 



3. Be — > P{y) transition form factors 

Be P{V) situations may be discussed in parallel by a corresponding displacement 
of inputs. We proceed to cope with the semileptonic transitions induced, respec- 
tively, by & ^ u, c and c — > d, s. However, since the related distribution amplitudes 
are not available from a QCD calculation, we give a model description for them, 
which is based on the harmonic oscillator potential ^1^. 

Our calculations, as in the B meson case, are limited to the regions where the 
OPE goes effectively, namely, < < ml — 2TOf,AQCD — 15 GeV^for 6-quark 
decays and < < — 2mcAQCD — 0.4 GeV^ for c-quark decays. The numerical 
results for the form factors at = are collected in Tab.l. It turns out that the 
calculated form factors can be parameterized excellently as 
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Table 1. The values of the form factors at ij^ = in comparison with the 
estimates in the three points sum rule (3PSR) (with the Coloumb corrections 
included) and in the quark model (QM). 



Mode 




/+(0) 


/-(O) 


MO) 


A+{0) 


A-(0) 


y(o) 


Be ^ 5c[lS] 


This work 
3PSr[13 

qmCHI 


0.87 
0.66 
0.76 


-0.87 
-0.36 
-0.38 


0.75 
0.63 
0.68 


1.69 
0.69 
0.66 


-1.69 
-1.13 
-1.13 


1.69 
1.03 
0.96 


Be ^ Bi*'' 


This work 

spsrCH 
qmCHI 


1.02 
1.3 
-0.61 


-1.02 
-5.8 
1.83 


1.01 
0.69 
-0.33 


9.04 
-2.34 
0.40 


-9.04 
-21.1 
10.4 


9.04 
12.9 
3.25 


Be ^ 


This work 

spsrCZI 
qmCHI 


0.90 
1.27 
-0.58 


-0.90 
-7.3 
2.14 


0.90 
0.84 
-0.27 


7.9 

-4.06 
0.60 


-7.9 
-29.0 
10.8 


7.9 
15.7 
3.27 


Be ^ D(*) 


This work 

3psr[I3 
qmCHI 


0.35 
0.32 
0.69 


-0.35 
-0.34 
-0.64 


0.32 
0.43 
0.56 


0.57 
0.51 
0.64 


-0.57 
-0.83 
-1.17 


0.57 
1.66 
0.98 



Table 2. Branching ratios (in %) of simileptonic Be decays into ground state charmo- 
nium states, and into ground charm and bottom meson states, in comparison with the 
result of 3PSR, QM, and the approach of the Bethe-Salpetcr equation. For the lifetime 
of the Be we take T{Be) = 0.45ps. 



Mode r/eCi^rieri^J/ipei^J/tpTi^ Deu\ Dtu\ D*eu\ D*tu \ Bei^ B* ei^ BsCi^ B*eu 



This work 1.64 


0.49 


2.37 


0.65 


0.020 


0.015 


0.035 


0.020 


0.21 


0.32 


3.03 


4.63 


3PSR 0.75 


0.23 


1.9 


0.48 


0.004 


0.002 


0.018 


0.008 


0.34 


0.58 


4.03 


5.06 


QM 0.98 


0.27 


2.30 


0.59 


0.018 


0.0094 


0.034 


0.019 


0.15 


0.16 


2.00 


2.6 


BSEli^ 0.97 




2.30 




0.006 




0.018 




0.16 


0.23 


1.82 


3.01 



Extrapolating the LCSR results with this parametrization, we can estimate the 
branching ratios of the simileptonic Be decays. The results are shown in Tab. 2 
together with those of other approaches, where we have input the following CKM- 
matrix elements: Vcb = 0.0413, Kb = 0.0037, Vcs = 0.974, Vcd = 0.224. For the 
6-quark decay modes, our results for the branching ratios are much larger than the 
corresponding those of 3PSR's. In these decays, the kinematical region is rather large 
and the main contributions to branching ratios should come from QCD dynamics at 
smaller recoil. The numerical values of the form factors always increase much faster 
with in the LCSR approach than in the simple pole approximation required in 
the 3PSR analysis, which accounts for the numerical discrepancy between the two 
approaches. While in the c-quark decays, where the kinematical region is narrow 
enough, the results of the two approaches are comparable to each other. 

4. summary 

In this paper we derive LCSR's for the form factors for semileptonic B{Bs,Bc) 
decays of current interest, which depend mainly on the leading twist distribution 
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amplitudes of the producing mesons by using suitable chiral currents. 

For better understanding the behavior of these form factors in the whole kine- 
matical range, a combined use is necessary of several different approaches which 
could be complementary to each other, for different kinematical regions. For exam- 
ple, in the B ^ D case one can use three different methods (LCSR, pQCD and 
HQET) for different kinematical regions to understand it in the whole kinematical 
range. 

Based on the models with the harmonic oscillator potential for the light-cone 
wave functions, we calculate the form factors for several important semileptonic 
Be decay modes in the their kinematical regions allowed by LCSR calculation. 
Extrapolating the LCSR results to the large regions, we derive the decay widths 
and branching ratios. 
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